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Abstract. Few body methods are used in many particle physics to describe 
correlations, bound states, and reactions in strongly correlated quantum sys- 
tems. Although this has already been recognized earlier, rigorous attempts to 
treat three-body collisions have only been done recently. In this talk I shall give 
examples and areas where few-body methods have been and might be of use in 
the future. 



1 Introduction 

Describing an ensemble of many particles (fcrmions/bosons) becomes challeng- 
ing and interesting as soon as interactions (e.g. Coulomb or strong interaction) 
are considered. Examples for Coulombic systems are ionic plasmas as they oc- 
cur in the sun and stars, liquid metals and electron-hole plasmas. Nuclear matter 
and the quark gluon plasma are examples of strongly interacting systems. Be- 
cause of the interaction it is not possible to treat even single particle dynamics 
without regarding effects of the other particles. Furtheron the system may be 
in equilibrium or out of equilibrium, depending on the boundary conditions 
imposed. 

The density temperature planes of matter and nuclear matter are shown 
in Figs, [l] and Q The phase diagram of nuclear matter turns out to be very 
rich. In particular the superfluid phase reflecting strong pairing is relevant for 
the structure of neutron stars M. At lower densities bound states occur. This 
part of the phase diagram may be accessed in the laboratory through heavy ion 
collisions at intermediate energies. The conditions for the formation of bound 
states are reached in particular during the final stage, where the nuclear density 
drops below the Mott density. 
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Figure 1. Density temperature phase 
diagram of matter. 
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Figure 2. Density temperature plane of 
nuclear matter. 



2 Theory 

Quantum statistics provides powerful methods to tackle the many particle sys- 
tems. Here I follow the Green function formalism ]2), which is convenient to 
introduce few-body methods. Let the Hamiltonian of the system be given by 

H(t)= ^ffoOOVl (t)Mt) +5^2(12, l'2')i>l(t)4(t)Tp 2 ,(t)Tp v (t), (1) 

1 12 

where ipi(t) denotes the Heisenberg operator of the particle with quantum 
numbers Si, fci, etc. for spin, momentum etc. The one particle Green function 
is defined by 



#1(1,1') = (TMWUt')) = Tr{p 2>i(i)Vp (t')}, 



(2) 



where averaging is due to the density operator pr,. For an open system in 
thcrmodynamical equilibrium the extremum condition for the entropy leads to 
the following expression for the quantum grand canonical density operator, 



Po = 



e -/3(H-nN) 
Tx{e-^ H -^ N )}' 



(3) 



The temperature (1//3 = T) and the chemical potential /1 are the correspond- 
ing Lagrange parameters. Using the Heisenberg equation for ipi results in the 
following equation for Qi 0] 

C?i(l, 1') - g[°\l, 1') - ]T g[°\l, l)iV 2 {12, 12) 02(12, l'2 + )| ti=t2 . (4) 

122 

The argument 2 + means that t 2 + = t 2 +0 + and 

(id tl -H (l))g[°\l,l')=Sn>. (5) 

Eq. |] shows already the basic problem of many particle physics, the hierarchy. 
To find a useful truncation of the n + 1 particle Green function from the n 
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particle one, some notion of the system is needed. Usually the hierarchy is 
truncated at the two particle level, assuming binary collisions only. Three- 
particle collisions have been treated at most approximately using Born (for 
Coulombic systems) or impulse approximation (for nuclear matter). This may 
not be sufficient, in particular, if explicit three-particle processes are considered 
(e.g. such as cluster formation, where a third particle is needed to achieve 
momentum conservation, etc.). 

Eq. |] may be formally decoupled by introducing the self energy 27(1, 1). 

£ 27(1, 1) Ox(l, 1') = - g iV 2 (12, 12) 02(12, 1'2+). (6) 

1 122 

In the simplest case the self energy may then be treated in mean field 
(e.g. Hartree-Fock) approximation, i.e. Q 2 — > — <?i(l,l') Qi(2,2') — 

Qi (1,2') Qi (2,1') viz. no two particle correlations (leading to an ideal gas of 
quasi particles). Using U 2 (12,12) = -U 2 (12,21) the self energy S HF (1,1') is 
given by 

Z HF (1, l')=iJ2 ^2(12, l'2)ft(2, 2+), (7) 

22 

which reduces to the standard expression, if a static potential is used, 
U 2 (12,1'2') = 5iv5 2 2>V 2 (12), since ft (2, 2+) = (V4V>2> = h is the one par- 
ticle distribution function for fermions. For a given potential Eq. [7] constitutes 
a self consistent problem to determine \x and /3, solved by iteration 




Figure 3. The nucleon self energy (left side). Devi- Figure 4. One particle distri- 
ation of the effective mass approximation from the bution function as a function of 
exact model result for small densities (right side), momentum k no = 17 fm -3 

As a further simplification useful in later applications effective masses may 
be introduced. The self energy for the nucleon in a nuclear medium of T = 
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10 MeV is shown in Fig. g along with its effective mass approximation e* = 
k 2 /(2m*) + U HF (k = 0). The distribution function f(k) is shown in Fig. £§ 
Calculations are done using a separable Yamguchi type potential (that will 
later be used for the three body calculations). Finally, analytic continuation 
of the Green function defined in Eq. leads to the Kubo-Martin-Schwinger 
boundary condition as e %tH — , i.e. 

Si(l,lO| tl =o = -e^Si(l,l')k=-i/3- (9) 

As a consequence the time like component of the Fourier transform is restricted 
to certain values only (Matsubara frequencies), Q^ 1 (l — > Gi(z„), where z v = 
mv / (3 + (i and v = ±1, ±3, . . . for fcrmions. 



3 Correlations 



A better treatment that goes beyond the quasi particle approximation is pro- 
vided, e.g. through the cluster approximation that include correlations ||. As 
a consequence the one particle spectral function Ai(u>), defined through 

J 2ir z v - lo 

which is given by A\{uj) = 2ttS(lo — e) for the quasi particle approximation, is 
more complicated, viz. 

A ( \ — 2InxE'(a;) . . 

l[UJ ' ~ \lo - E - + [lmZ(Lu)} 2 ' ( ' 




Figure 6. Temperature density 
Figure 5. The spectral function Ai(u>) for dif- plane for nuclear matter as a func _ 
ferent temperatures and densities including two tion of the total density . The amount 
particle correlations. The deuteron bound state Q f corre i a tes density is given in per 
is recognized as the left shoulder. cent of the total density . 



As an example Figure [| shows the spectral function using the full two 
body t-matrix to describe the correlations [H . As a consequence nuclear matter 



5 



cannot be considered as a system of independent quasi particles but for a large 
part it is correlated up to full pairing in the superfluid phase. This is depicted in 
Fig. U where the dashed lines show equal contribution of the correlated density 
to the total density ||. The basis to treat correlated densities is provided 
by a generalization of the Beth-Uhlenbeck approach ||. The nuclear density 
n = n((i, T) is given by 



7lf r , 



2n 2 + 3n 3 + 



(12) 



where n 2 ^ denotes the two, three-particle correlations, present as 
bound/scattering states. In first iteration these correlations may be treated 
on the basis of residual interactions between the quasi particles. The exact two 
particle equations to be solved are known as Bethe-Goldstone or Feynman- 
Galitski equations depending on some details. In ladder approximation the 
equation for the two body Green function reads 



G 2 (z) 



/1/2 - /1/2 , /1/2 — /1/2 



£2 



£1 - £2 



V 2 G 2 (z), 



(13) 



where / = 1 — /. Introducing the two-body i-matrix in a standard fashion both 
bound and scattering states have been solved 
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Figure 7. The deuteron binding energy 
as a function of the nuclear density for 
T = 10 MeV. P denotes the relative mo- 
mentum between the deuteron and the 
medium. 



Figure 8. The triton binding energy as 
a function of nuclear density for different 
temperatures T. The triton rests in the 
medium. The respective Nd continua are 
also indicated. 



As an example the deuteron energy is shown in Fig. 0. The Mott density 
is defined through the condition E c i = 0. Note that the Fermi functions fi for 
particle 1 etc. depends on the relative momentum P between the deuteron and 
the medium. The respective three-particle Faddeev type equation for the Green 
function has been given in Ref. Q| , 



G 3 (z) 
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Eq. [l4| has been rewritten using the AGS approach M . The resulting in- medium 
AGS equations are 

U a/ 3 = {l-5 aP ) — — +2,^'— ; —U 7fl . (15) 

\z — £1 — £ 2 — £3/ ~r z — £1 — £2 — £3 

where Tg 7 " 1 is the solution of the in-medium two-body problem, e.g. for 7 = 3 

T 3 (3) = (1 - f 3 + g(e 1 + e 2 ))- 1 V2 + V 2 hh ~ hh T 3 (3) . (16) 

z — £1 — £2 — £3 

and the Bose function is given by g(ui) = 1/ '(e^ -2 ^ — 1). This equation has 
been solved for the Nd reaction relevant for deuteron formation ^| (sec 
below), assuming f(e) 2 <C /(e), compare Fig. ||. Recently, the triton bound 
state equation has also been solved |l(J. The triton binding energy is shown 
in Fig. ||. The deuteron as well as the triton binding energies weaken if the 
nuclear density is increased until the Mott density is reached. This tendency is 
dominated by the Pauli blocking of the surrounding medium. 

The four-nucleon correlation is believed to play a significant role for lower 
densities and temperatures. Exploratory calculations using a simple variational 
ansatz for the 4 He wave function predict an a condensate/quartetting on top 
of the deuteron condensate/triplet pairing that leads to superfluidity p] ]. 

Bose systems behave quite differently with respect to the occurrence of 
bound states. The bose functions enhance the effective residual interaction 
that might lead to an "opposite Mott effect" , i.e. existence of bound states and 
also pairing (condensate) even if no bound state exist for the isolated case. An 
example is provided by a pion gas, were a pionic condensate may occur fl^l 
that are discussed, e.g. in the context of neutron stars jjl. 



4 Reactions 



Nuclear reaction rates play an essential role in the formation of stars like the 
sun. The standard solar model is based on binary collisions. Recently, triple 
reactions, e.g. e + 3 He + a S Be + e to be compared to 3 He 4- a — > 7 Be + 
7 have been investigated and found to be rather small in plasmas at solar 
conditions |r|. However, note that triple collisions are mostly non-radiative 
and that they may be more important for other stars than the sun or at the 
early universe |i3fl . 

Another example are dense ionic plasmas, where the ionisation rate depends 
on three-particle reactions that are presently treated in Born approximation. 
Since the residual interaction is Coulombic this may be considered a good 
approximation and it reproduces the experimental results for hydrogen like 
plasmas |Q. For higher ionized plasmas this might not be the case and the 
application of Faddeev like methods may be in order. These will be sketched 
in the following for nuclear matter. 
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The generalized quantum kinetic Boltzmann equations for the nucleon 
/n(j>, t) (momentum p) and deuteron f d (P, t) (momentum P) distribution func- 
tions |u| 

f N (p,t) = (a' Np a Np ) = Tr{p(t)a' Np a Np } 
fd(P,t) = (bl P b dP ) (17) 
are coupled and read 

/iv(p,t) = -T> N (p,t) +X N (p,t) 

f d (P,t) = -V d (P,i)+I d (P,t), (18) 

The first term reflects the so called Vlasov term and is related to the mean 
field. The second term is the collision term that is responsible for equilibration 
of the system. The explicit form of the integral In(p, t) is 

l N {p,t) =!>{p,t)f N {p,t) -I<{p,t)f N (p,t), (19) 

where, e.g. T^(p,t) is given by 

ZjvG 5 )*) = [ dkd k 1 dk 2 \(kp\T NN _ NN \k 1 k 2 )\ 2 f N (k 1 ,t)f N (k 2 ,t)f N (k,t) 

dkdki dk 2 dk 3 \(kp\U Nd^NNN | h k 2 k 3 ) | 2 

xf N (k u t)f N {k 2 , t)f N (k 3 , t)f d (k, t) + ... (20) 

The solution of this equation are the distribution functions /at and f dl that how- 
ever also appear in the transition matrices Tnn->nn, U Nd^NNN, etc. Therefore 
a full solution of this equation is a difficult problem. For small fluctuations from 
the equilibrium distributions the equations may be linearized in the framework 
of linear response theory. The binary collision approximation may the transition 
matrix elements depend on the equilibrium distribution only and the results 
of the previous sections can be utilized. Even this in-medium dependence has 
hardly been considered in modeling of heavy ion collision fl?| . Also, for three- 
nucleon collisions so far only the impulse approximation has been used [[L7| . 
Here we solve the in-medium Faddeev type equation that includes the Hartree- 
Fock self energy shift and the Pauli blocking in a consistent way, Eq. [l4| The 
resulting break-up cross section for a typical temperature of T = 10 MeV and 
densities below the deuteron Mott density is shown in Fig. ^[ For a comparison 
of the quality of the model the isolated cross section along with the experi- 
mental data are shown in Fig. |l(| From inspection of Fig. ^ we see that 
the in-medium cross section is significantly enhanced compared to the isolated 
on. The threshold is shifted to smaller energies, which is because the binding 
energy of the deuteron becomes smaller. We observe that for higher energies 
the medium dependence of the cross section becomes much weaker. 

Though the change of the cross section looks dramatic, the quantum Boltz- 
mann equation still has some additional medium dependence that may change 
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Figure 9. In-medium break-up cross 
section at T — 10 MeV. Isolated cross 
section is shown as solid line. Other lines 
are due to different nuclear densities. 



Figure 10. A comparison of the to- 
tal, elastic, and break-up cross sections 
nd — > nd, nd —> nnp with the experi- 
mental data of Ref. ^[ . 



this effect in the observables /n and fd- Within the linear response theory it 
is possible to calculate the chemical reaction time due to the break-up process. 
For small fluctuations 5f(t) linear response leads to 

d t Sf d e&ction (P, t) = —Sf d caction (P, t) (21) 

TP 

where the "life time" of deuteron fluctuations has been introduced, 

_ 1 4 f 2 

t a = y I dk N dk 1 dk 2 dk 3 \{kp\U Nd ^ NNN \kik2k 3 )\ /1/2/3/s 2irS(E - E ). 

(22) 

which can be related to the break-up cross section given in Eq. |j| For low 
densities the life time (as a function of the deuteron momentum P) and the 
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Figure 11. Fluctuation time 
of the deuteron distribution as a 
function of the deuteron momen- 
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inverse life time, i.e. the width, at P = along with the deuteron binding energy 
for comparison is shown in Fig. [H] jj| . These times have to be compared to the 
approximate duration of the heavy ion collision of about 200 fm. 

5 Conclusion 

The treatment of correlations and triple collisions in non-ideal many particle 
quantum systems opens up a new field for few-body methods. The examples 
shown have been mostly from nuclear physics. However, applications are pos- 
sible for stellar plasmas to improve the description of the basic quantity, which 
is the spectral function and to include correlations into the equation of state. 
In the laboratory the ionisation rate of dense ionic plasmas is determined by 
three-particle collisions. A description is presently restricted to hydrogen like 
plasmas, where the Born approximation for the three-particle reaction is suffi- 
cient. Typical applications in nuclear physics are related to heavy ion collisions, 
here in particular the formation of light clusters such as deuterons, helium-3, 
tritium, and alpha particles. The conditions are satisfied during the final stage 
of the heavy ion collision, where a temperature of T ~ 5 ... 10 MeV may be 
meaningful and the densities are below the Mott densities of cluster formation. 
The results are also relevant for the equation of states of neutron stars. 

The approach given here follows the quantum statistical description as it 
provides a rigorous, systematic treatment of many particle systems. The major 
approximation utilized is the cluster expansion to decouple the infinite hier- 
archy of equations (Green functions or kinetic equations). This approximation 
clearly goes beyond the quasi particle picture as it includes the residual inter- 
actions in a systematic way. This is done rigorously using few-body methods. 
To this end few-body equations have to be substantially generalized. Presently, 
these equations resemble an RPA structure however extended to finite 
temperatures. 

The validity of this already ambitious approach has to be checked by facing 
experimental results. This seems easier for ionic plasmas, e.g. to calculate the 
ionisation rate, or for electron-hole plasmas in the context of exciton formation. 
Testing the validity of the approach for nuclear physics needs a handle of heavy 
ion collisions. Presently one relays on numerical simulations of the complicated 
dynamics of a heavy ion collision, which is subject to discussions by its own. 
Typical heavy ion simulation codes that require microscopic input are based 
on e.g. a Boltzmann-Uehling-Uhlenbeck treatment or on quantum molecular 
dynamics. The results presented here may however also be relevant for standard 
nuclear physics, e.g. electron scattering off heavy nuclei, when correlations are 
considered, and one therefore needs to go beyond the quasi particle picture. 

Acknowledgement. It is a pleasure to thank my colleagues in Rostock who 
provided me with some material presented during the talk. 
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